rm(list = ls())
load("dataBJPOLS_weekdays.RData")

cases.weekday <- last.cases
round(c(mean(cases.weekday$pti, na.rm = T), sd(cases.weekday$pti, na.rm = T)), 2)
round(c(mean(cases.weekday$daysjail_t > 365, na.rm = T), sd(cases.weekday$daysjail_t > 365, na.rm = T)), 2)
round(c(mean(cases.weekday$daysjail_t, na.rm = T), sd(cases.weekday$daysjail_t, na.rm = T)), 2)
round(c(mean(cases.weekday$age, na.rm = T), sd(cases.weekday$age, na.rm = T)), 2)
round(c(mean(cases.weekday$female, na.rm = T), sd(cases.weekday$female, na.rm = T)), 2)
round(c(mean(cases.weekday$race4 == "B", na.rm = T), sd(cases.weekday$race4 == "B", na.rm = T)), 2)
round(c(mean(cases.weekday$race4 == "W", na.rm = T), sd(cases.weekday$race4 == "W", na.rm = T)), 2)
round(c(mean(cases.weekday$race4 == "H", na.rm = T), sd(cases.weekday$race4 == "H", na.rm = T)), 2)
round(c(mean(cases.weekday$any_drug, na.rm = T), sd(cases.weekday$any_drug, na.rm = T)), 2)
round(c(mean(cases.weekday$any_weapon, na.rm = T), sd(cases.weekday$any_weapon, na.rm = T)), 2)
round(c(mean(cases.weekday$any_armed_rob, na.rm = T), sd(cases.weekday$any_armed_rob, na.rm = T)), 2)
round(c(mean(cases.weekday$any_prior_case, na.rm = T), sd(cases.weekday$any_prior_case, na.rm = T)), 2)
round(c(mean(cases.weekday$bond_amt_n_max, na.rm = T), sd(cases.weekday$bond_amt_n_max, na.rm = T)), 2)
